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ABSTRACT 

The European VLBI Network (EVN) has been used at two epochs in 2003 and 2009 to obtain 
multi-frequency high-resolution images of the merger galaxy NGC 6240 in order to study the 
radio properties of all compact high-brightness components in the galaxy. Our observations at 
milli-arcsecond resolution detected the northern and southern nuclei and two radio components, 
which we interpret as long-lived luminous supernovae associated with the circum-nuclear star- 
burst activity at the southern nucleus. The new VLBI data support the presence of an active 
galactic nucleus (AGN) together with starburst activity at the southern nucleus and provides 
some evidence for an AGN at the northern nucleus. The two nuclei both display an inverted 
spectrum at lower GHz frequencies. The spectrum of the southern nucleus indicates thermal 
free-free absorption on parsec scales, consistent with the presence of an AGN. 



Subject headings: galaxies: active — galaxies: individual (NGC 6240) 
continuum: galaxies — supernovae: general 



Jaxies: starburst 



radio 



1. Introduction 

The merging galaxy NGC 6240 has a large far- 
infr ared (FIR) luminos ity of Lfir = 3.5 x 10 11 
L ( Yun fc Carilli|[2002T) and belongs to the cate- 
gory of luminous inf rared galaxies (LIRGs; Lt r = 
10"-10 12 L ) (see ISanders fc Mirabell fl996l . for 
a review). The large luminosity of these LIRGs 



and Ultra-Luminous Infrared Galaxies (ULIRGs; 
L m = 10 12 -10 13 L©) is likely to be domi- 
nated by starburst activity induced by galaxy- 
galaxy interactions, where the UV radiation from 
massive star formation is reprocessed to far- 
infrared radia ti on in t he dusty environmen t (e.g. 
ISanders et al] Il988i: ISanders fc Mirabell Il99 
19971) . The NGC 6240 merging sys- 
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tern has been descri bed at optical wavelengths 
(|Fried fc Schuldll983l) and is thought to host an 
active galactic nucleus (AGN) in its nuclea r regio n 



(e.g., iDepov. Becklin. fc Wvnn- Williams! [l986) 



The two galactic nuclei found in the central re- 
gion have a projected separation between 1".5 



1 



and 1".8 (0. 714-0.856 kpc) at X - ray and radio 
wave l engths (Fried fe Schulall983t iBeswick et al 

iGalli 



2001; Komossa ct 



all 1200? 



2004[ iMax et al.l (2007). 



free absorption (FF A) or possibly synchrot r on-se lf 



iGallimore fc Beswick 
Because the north- 
ern nucleus, with the highest absorbing column 
densities, lies well behind the s outhern nucleus 
( Baan. Haeiwara fc Hofnerll2007h . ;he separation 
varies as the intervening hot (obscuring) dust 
becomes optically thinner at longer wavelengths 



(jMax et al.l 120071) . X-ray observations at 2—10 



keV reveal the presence of iron line emission at 
6.4 keV at both nuclei, which is most prominent 
in the southern nucleus and classifies this as a 



binary AGN ([Komossa et al.l 120031 ). 

Thermal CO emission reveals a significant mass 
conc entration that is centered between the two nu- 
clei dTacconi et al.lll999t iBrvant fe Scovillej [l999t 
iNakanishi et al.ll2005Hlono et al|2007tl andlargel y 



consists of interstellar gas (IScoville et all l2000l) 



This intervening molecular gas concentration re- 
sults from a superposition of disk gas from the two 
merging galaxies, which is confirmed by the struc- 
ture and peak location of the OH absorption struc- 
ture against t he extended radio emission of t he 
nuclear region ([Baan. Haeiwara fe Hofnerll2007l) . 

Radio observations at centimeter wavelengths 
of the central 2" region of NGC 6240 using the 
NRAO Very Large Array display the two nuclei 
and an extended northern structure (N3) that 
is part of the large-scale lo op structure to the 
west of the nuclear sourc e ( Colbert et al. 1 119941: 
Baan. Haeiwara fc Hofnerl 120071 ). A large scale 



diffuse component envelops these three structural 
components. The nuclear radio components re- 
main unresolved at higher resolutio n at 5 GHz 



using MERLIN (|Beswick et al.ll200ll ). The loop 



structures including N3 have been resolved away 
at this resolution and the southern nucleus, Nl 
shows east-west structure at its lowest contours 
(PA = 60°) and a northern extension (PA = -30°). 

At milli-arcsecond resolution using the NRAO 
Very Long Baseline Array (VLB A), the northern 
nucleus N2 also displays a multi-component east- 
west stnicture_£i ; rg^estmga core with a two-sided 



jet (jGallimore fc Beswick I I2004T) . The southern 



nucleus has a peak central structure surrounded 
by weak additional components possibly related 
to starformation. Both high-brightness nuclei 
show an inverted spectrum at lower frequency 
with a spectral turn-over or flatness due to free- 



absorption (SSA) (jGallimore fc Beswick 1 120041) 



Strong H 2 Q maser em ission was discovered (e.g., 
Haeiwara et al.l [2002 ) at the exact center of the 
south ern nucleus ([Haeiwara et al.ll2003l ; lHagiwara 
2010), which may be related to a dense molecular 
gas around a central engine. 

Here, we present European VLBI Network 
(EVN) observations with increased surface bright- 
ness sensitivity to image the compact sources at 
increasing frequencies. New EVN observations at 
1.6, 2.2, 5.0 and 8.4 GHz were made to study the 
source structure, the overall spectral energy dis- 
tribution, and to check for consistency with the 
earlier VLBA results. In Section 2 we present the 
new observations, the data analysis, and the obser- 
vational results, including the identification and 
morphology of all detected compact radio sources. 
In Section 3 we discuss the nature of these sources 
and compare NGC 6240 with other active galaxies. 
Finally, in Section 4 we summarize the discussion. 

Cosmological parameters of H = 73 km 
s~ 1 Mpc~ 1 , = 0.73, and = 0.27 are adopted 
throughout this article, which results in a lumi- 
nosity distance of NGC 6240 of 103 Mpc at z 
= 0.0245 and an angular conversion of 1 arcsec 
equaling 476 pc in linear scale. The designa- 
tion of radio compo nents follows that defined in 
Colbert et al.l (| 19941 ) with Nl being the southwest- 
ern nucleus and N2 being the northeastern nuclei. 
T he two nuclei were called resp ectively S and Nl 
"~ (1200 
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2. Imaging the double nuclei in NGC 6240 
2.1. Observations and imaging 

The European Very Long Baseline Interferom- 
etry Network (EVN) has been used to observe the 
nuclear region of NGC 6240 in order to investigate 
the nature of the active galactic nuclei and nuclear 
starbursts in the galaxy. The stations used for 
these observations at two epochs in 2003 and 2009 
were Effelsberg (Ef), Hartebeesthoek (Hh), the 76- 
m Lovell telescope (Jb) at Jodrell Bank, Medic- 
ina (Mc), Noto (Nt), Onsala (On), Shanghai (Sh), 
Torun (Tr), Urumqi (Ur), and the Westerbork 
(Wb) phased-array. A summary of the observa- 
tional details is presented in Table [T] All observa- 
tions were conducted in phase-referencing mode at 
all frequency bands using the nearest calibration 
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source J1651+0129 (centered at R.A.(J2000) = 
16 h 51 m 03 s .6620, Dec.(J2000) = 01°29'23".458) at 
1.01° away from NGC 6240. Only the Lovell tele- 
scope did not participate in the phase-referencing 
part of the observations due to its lower slewing 
rate. 

The first epoch observations at 1.6 GHz and 5.0 
GHz were made on 30 October and 10 November, 
2003 and utilized a 32 MHz total bandwidth with 
2 bit sampling in dual-circular polarization mode 
for a total data recording rate of 256 Mbit s _1 . 
Observations of 2.5—3 minutes on the target were 
interleaved with observations of 1.5—2 minutes on 
the phase-referencing source, which gave about 4.2 
hours total on-source time. 

Second epoch observations at 5.0 GHz, 2.2 
GHz, and 8.4 GHz were made on 16, 17 June, 
2009. The observations at 5.0 GHz were made in 
dual-circular polarization mode utilizing 64 MHz 
total bandwidth with 2 bit sampling for a to- 
tal data recording rate of 512 Mbit s _1 . Three 
minute scans on the target were alternated with 
one minute scans on the phase-referencing source, 
which resulted in 5.7 hours total on-source time. 
The observations at 2.2 and 8.4 GHz were made 
in single polarization mode with the S/X dual- 
frequency receiving system with a bandwidth of 
128 MHz at each frequency and 2 bit sampling 
for a total recording rate of 512 Mbit s _1 . Three 
minute scans on the target and one minute scans 
on the phase-referencing source were alternated 
during the observation, which resulted in 2.9 hrs 
total on-source time at each frequency. All ob- 
servations were done with the same reference po- 
sition near the southern nucleus of NGC 6240 
(R.A.(J2000) = 16 h 52 m 58 s .8903, Dec.(J2000) = 
02°24'03".339). 

The data were correlated using the EVN MklV 
correlator at the Joint Institute for VLBI in Eu- 
rope (JIVE) and the output visibility data were 
averaged to 1-2 seconds. Sixteen frequency chan- 
nels of 500 kHz were used across each 8 MHz band. 
Data calibration and imaging were carried out in 
the NRAO Astronomical Image Processing Sys- 
tem (AIPS). After initial data editing, the delay 
error corrections resulting from ionospheric effects 
were applied to all visibility data. A-priori am- 
plitude calibration was performed using the sys- 
tem temperature and gain information provided 
for each telescope. The delays and fringe rates 



for the calibration sources 3C 273 or 3C 345, and 
J1651+0129 were applied to the target source, 
NGC 6240. Using the measured flux densities 
of the phase-referencing source J1651+0129 at 2 
epochs as a basis, the relative accuracy of the tar- 
get flux densities is estimated to be ~10%. This 
error estimate includes the effects of bandwidth 
and time smearing and ionospheric delays, which 
are estimated to be several percent for these EVN 
observations. 

Cleaned images of 512 x 512 pixels were pro- 
duced, one centered near the southern nucleus and 
one at the intensity peak position of the north- 
ern nucleus offset by 0".5051 in right ascension 
and 1".4279 in declination from the phase cen- 
ter. The noise levels of 1—2 times theoretical 
noise and synthesized beams from CLEAN-ed im- 
ages are summarized in Table O Self-calibration 
could not be applied successfully because of the 
low signal-to-noise ratio of the images. Final im- 
ages have been phase-referenced to the nearby cal- 
ibrator J1651+0129. Because of residual errors in 
the calibration of the phase-reference sources, the 
associated structure at the lowest contours (< 3 a 
or 4 a) in the source images may not be reliable. 

The positions of the radio components and their 
errors at the peak brightness at each frequency and 
both epochs are listed inTable El These compos- 
ite position errors have been estimated from the 
theoretical t hermal noise errors of interferometer 
phase (e.g.. iThompson et alJ ll986). the system- 



atic errors of the phase-referencing VLBI observa- 
tions (cphr), and the errors of the absolute position 
of the phase-referencing source (cr rc f), using the 

quadrature equation: 



Here 6 represents the synthesized beam size and 
SNR represents the signal-to-noise ratio. The val- 
ues of CTphr (Aa cos<5 = 0.060 mas, AS = 0.020 mas) 
are obtained from the s imulated astromet ric accu- 
racy of the EVN array (jPradel et al.ll2006l ). adopt- 
ing a separation between calibrator and target 
of 1.0° at position angle 152.2° and at an ap- 
proximate source declination of 0°. The value of 
OY B f is 0.58 (mas) fro m the VLBA Calibrator List 
(jKovalev et alj|2007l ). Table H lists the flux den- 
sities of the radio components at each observing 
frequency, which also includes two weaker compo- 
nents appearing in the earli er VLBA observations 
(see Figures 2, 3, and 4 in iGallimore fc Beswick 
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20041 ) . Source sizes of these components in Table [5] 
are measured from the hnal images using the AIPS 
task JMFIT. Assuming Gaussian distributions for 
sources, the upper limits of the source sizes have 
estimated accuracies of 3—30% depending on ob- 
serving frequency and epoch. Table [6] presents the 
properties of the radio sources detected in our ob- 
servations. 

2.2. Source identification and morphology 

The new VLBI maps obtained from 1.6—8.4 
GHz display four compact high-brightness emis- 
sion components as identified in Figure [TJ the nu- 
clei, Nl (south) and N2 (north ) identified in earlier 



intcrferometric observations (Beswick et al. 2001 



Hagiwara ct al. 2003) and two new components 
RSI and RS2 near the southern nucleus. The ra- 
dio source RSI was first identified southwest of 
the southern nucleus at 5.0 GHz in October, 2003 
and it is detected at 5.0 and 8.4 GHz in June, 
2009 (see Figure [2]). The other radio source RS2 
was identified northeast of the southern nucleus 
at Epoch 1 and was more prominent at 1.6 GHz 
than at 5.0 GHz (see Figure 0]). At Epoch 2 
the source RS2 was clearly detected at 5.0 GHz 
but not at 2.2 GHz with an upper limit flux 
density of 2.75 mJy beam -1 . All diffuse emis- 
sion in the central 2" region, obse rved in earlier 
Very Large Array observations (e.g.lCoibert et al 



1994 ; Baan. Hagiwara fe Hofner 20071) . has been 
resolved out. 

The VLBA observation at 8.4 GHz of August 
1999 did not detect the tw o nuclei at a 5a upper 
limit of 0.65 mJy beam -1 ([Gallimore fc Beswick 



2004), but the northern nucleus was detected at 
8.4 GHz in June 2009 with similar sensitivity, 
which suggests some intensity variability of the 
northern nucleus. Both nuclei were detected at 
5 GHz at both Epochs. Their peak positions in 
the EVN data (Table [3]) agree at all frequencies 
with earlier VLBA measurements using estimated 
uncertainties of 0.6 mas in right ascension and 1. 5 
mas in declination (Galli more fc Beswick 2004). 
The 5.0 GHz maps in Figures [2] and [3] display the 
two nuclear components Nl and N2 obtained at 
Epoch 1 and Epoch 2, respectively. The southeast- 
ern nucleus, Nl was not detected at the highest 
frequency of 8.4 GHz with an upper limit of 0.55 
mJy beam - 1 (Figure [2|). Individual source maps 
indicate that the sources Nl and N2 are resolved 



but that the minor axis of RS2 northeast of Nl is 
unresolved (Figure S]). The peak of RSI is sepa- 
rated from Nl by 23.3 mas at 5.0 GHz at a posi- 
tion angle of ^35°, corresponding to 11.1 pc at the 
distance to NGC 6240. At the lowest frequency of 
1.6 GHz, Nl and RSI are not spatially separated. 
Also, RSI was not detected at 2.2 GHz at Epoch 2 
possibly because of the lower signal-to-noise ratio. 

The source RS2 is located northeast of the 
southern nucleus in our 1.6 GHz EVN map of 
Figure [4] It appears that the peak position 
of RS2 at 1.6 and 5.0 GHz lies close to the 
compact component SI at 0".3 northeast from 
Nl in the 1.7 and 2.4 GHz VLBA observations 
( Gallimore fc Beswick 20041 ). On the other hand, 
RS2 has an offset of « 10—15 mas at 1.6 GHz 
or 5.0 GHz relative to SI at 1.7 GHz or 2.4 GHz, 
which is not within the estimated positional errors 
(Table[3]). Therefore, RS2 is a different source to 
SI in the earlier VLBA data. We interpret both 
components to be part of the circum-nuclear star- 
forming region that extends over 10 mas or ~5 pc 
around Nl. 

The structure of the northeastern nucleus N2 
at 5.0 GHz (Figure [3]) shows a northwest ex- 
tension of 0".02 (w 9.5 pc) that is consis- 
tent with the earlier VLBA maps at 1.7 and 2.4 



GHz (jGallimore fc Beswick 1120041 ). Although this 



structure at lower contours may not be reliable be- 
cause of incomplete phase and amplitude calibra- 
tion, an east-west structure may be explained as a 
core-jet structure typically seen in high-luminosity 
Seyfert nuclei. The northern nucleus itself remains 
unresolved at both 1.6 and 2.2 GHz. Only at the 
highest frequency of 8.4 GHz, N2 has an east- 
west elongation that is consistent with the earlier 
VLBA data (Figure [3]). 

3. Discussion 

Our EVN data present the detection of the dou- 
ble nuclei at higher frequencies of 5.0 and 8.4 GHz, 
where earlier VLBA data failed to detect these nu- 
clei. These results combined with earlier VLBA 
studies, allow study of the nature of the nuclei 
at milliarcsecond-scale structure. While our study 
focuses on the physics of the two nuclei, the de- 
tection of two new components is a key to under- 
standing the activity of NGC 6240. 
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3.1. The northern nucleus N2 

The northern nucleus N2 shows an inverted 
spectrum rather than a steep spectrum that would 
be more typical for a radio jet (Figure [3]). The 
structure of N2 at 5.0 GHz is partially resolved 
and shows an east-west extension at its lower con- 
tours that still remains questionable because of 
the insufficient amplitude calibration (Figure [3]) . 
Changes in source structure may also result from 
a different (u,v) coverage for the two epochs. 
The detection of the X-ray emission towards N2 
(jKomossa et al.l 120031 ) is less convincing than the 
definite detection towards Nl. Nevertheless, the 
radio morphology of N2 resembles that of a core- 
jet st ructure seen in Seyfer t nuclei, such as NGC 
1068 (jGallimore et al.ll2004l ). but this requires fur- 
ther verification. By contrast with Arp 299, the 
VLA observations of the H 2 maser in NGC 6240 
have not detecte d any emission towards N2 dur- 
ing s even years ( Hagiwara et al. 20031: Hagiwara 
2010h . 



Baan. Hagiwara fc Hofnen 120071 ) and may origi- 
nate in the receding side of a co mpact rotating 



3.2. The southern nucleus Nl 

The southern nucleus Nl appears partially re- 
solved in an east- west direction (Figure[2p,2c) but 
its flux density at 5.0 GHz has not changed signif- 
icantly during 6 years. Contrary to the spectrum 
of N2 (see discussion below) , the spectrum of Nl 
with missing Epoch 2 detections at 2.2 and 8.4 
GHz may not be inverted (Figure [5]), although its 
spectral shape (including upper limits) strongly 
indicates a spectral turnover around 2.0 GHz. The 
non-detection of Nl at 8.4 GHz suggests a rela- 
tively flat AGN-like spectrum and a lower limit for 
the spectral index a < 0.55 between 5.0 and 8.4 
GHz, using oc v a with being the flux density 
at frequency v. In addition, the radio brightness 
temperature > 10 6 K (Table E]) and the detection 
of strong hard X-ray emission and neutral iron 
line emission would argue for the presence of an 
AGN in Nl. Considering the available evidence, 
the southern nucleus may host both an AGN and 
a (circum-)nuclear starburst region. 

The H2O maser features in NGC 6240 nearly 
coincide with the continuum peak of the souther n 
nucleus ( Hagiwara et al. 2003 ; Hagiwara 2010() . 
The narrow H2O maser lines are redshiftcd 
by ~ 200—300 km s _1 with respect to the 



mole cular disk at the nucleus (e.g., iMivoshi et al 
1995J), which would also support the presence of 
an AGN. Alternatively, these narrow lines dis- 
tributed over 120 km s _1 could be explained by 
wind maser emission as observed in the Circinus 



galaxy (jGreenhill et al.ll2003l ). 

The southern nucleus Nl of NGC 6240 has char- 
acteristics that are co mparable with tho s e of com- 
pone nt A of Arp 299 jNeff et all 12004 lUlvestad 
20091 ). The 8.4 GHz r adio power of Nl of (P N1 = 
3.7 x 10 22 W Hz" 1 ) (jColbert et al.lll994l) is only 
a factor of two more than the 8.4 GHz power of 
nucleus A of Arp 299 (P A = \-8 x 10 22 W Hz" 1 and 
P m =2.9 x 10 21 W Hz" 1 ) (|Neff et al.ll200l . The 
compact source A of Arp 299 has been resolved 
into discrete sources at milliarcsecond-scale by re- 
cent EVN observations and one of these sources 
(Al) is the AGN cand idate with a flat spectra l 



index of -0.13 ± 0.11 (jPerez-Torres et alj BoiO) 



systemic velocity of Nl (e.g., Hagiwara 2010; 



Five bright and compact radio sources lie within 
10 pc and a re candidates of Type-II young radio 
supernovae (jNeff et al.l 120041 ) . Four of these five 
radio sources have flat or inverted radio spectra 
between 2.2 and 8.4 GHz, very similar to RSI 
and RS2 near Nl in NGC 6240. The luminous 
H2O maser emission i n Arp 299 is found towards 
both nuclei A and Bl (jHenkel et al.ll2005l ) and at 
A it coinci des with its conti nuum peak at VLA 
resolution ( Tarchi et al. 2010l ). While this spatial 
coinci dence may also suggest an AGN association 
(e.g., iHagiwaral 120071 ) . the maser in Arp 299 is 
blueshiftcd relative to the systemic velocity and 
is likely associated with an outflow in a starform- 
ing environment. 

3.3. Spectral analysis of Nl and N2 

In order to understand the observed spectra 
of the southern and northern nuclei, model fit- 
ting may be used to explain the frequency turn- 
over at lower frequencies using a power-law plus 
pure free-free absorption ( FFA) or synchrotro n 
self-absorption (SSA) (e.g.. iKameno et aL 2000). 
First, the FFA model for a nuclear spectrum is 
described by: 

S v = S (v/1.0) ao cxp{-r ff (^/1.0)- 2 - 1 }, (1) 

where v is the frequency in GHz, So is the unob- 
scured synchrotron flux density in mJy, «o is the 
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optically thin non-thermal spectral index, and rg- 
is the opacity at 1.0 GHz. This model describes 
all FFA in the foreground to the synchrotron emis- 
sion source and assumes that all source compo- 
nents a£e_j3ubject tothe same foreground opacity 
(e.g JParra et alj2007t ). The second model of pure 
SSA is described as follows: 



S v 



exp{- 



5 }] 



(2) 



where v is the frequency in GHz, Sqt s will be close 
to the flux density, if the SSA coefficient t s << 1, 
and ao is the optically thin non-thermal spectral 
index. 

The parameters of the FFA and SSA spectral 
fits for Nl and N2 and the reduced \ 2 values (per 
degree of freedom) have been summarized in Ta- 
ble \7\ and the fitted curves are shown in Fig- 
ure [5] Both the FFA and SSA models produce 
reliable fits for Nl with low values (< 1.5) for 
the reduced % 2 with the FFA model being slightly 
worse but does not make good fits to N2. How- 
ever, since there is no compact source in Nl or N2 
having a high brightness temperature, any syn- 
chrotron self-absorption in these nu clei is unlikely 
( Kellermann fc Paulinv-Toth 1969). Using con- 
ventional energy equipartition arguments at lower 
radio frequencies, a turn-over frequency of 2.0 — 
3.0 GHz for a self-absorbed source requires an 
equ ipartition brightn ess temperature of ~10 11,95 



K (jReadhead 1 119941 ). The observed brightness 
temperatures for our sources are only ~ 10 6 — 10 7 
K. Better constraints of the nuclear properties of 
Nl in terms of model fitting requires more com- 
plete flux density measurements. 

While for N2 neither the pure FFA nor the 
pure SSA model can perfectly explain the spec- 
tral bending at 2—3 GHz (see Table [7] and Fig- 
ure E|) , the observed bending does indicate that 
free-free self-absorption by ionized foreground gas 
in a starburst environment is relevant. Since there 
is also no hint of self-absorption i n the MER- 
LIN spectrum (jBeswick et al. 2001 ). the size of 
such an absorbing medium must be very small 
and between about 9.5 and 25 pc in linear size 
(jGallimore &: Beswick 1120041) . Similarly, the pres- 
ence of an AGN in N2 cannot be confirmed from 
our EVN data, except that the core-jet like struc- 
ture, as also seen in the 1.7 and 2.4 GHz VLB A 
images, could support the presence of an active 
nucleus. 



3.4. Radio sources near the southern nu- 
cleus RSI, RS2 

The radio source RSI has a radio power at 8.4 
GHz of 4.1 x 10 21 W Hz -1 and a spectral index of 
"5.0-8.4 ~ 0.34 in our Epoch 2 observations. The 
radio power of RS2 at 5.0 GHz (ai. 6-5.0 ~ -0.91) 
is 7.5 x 10 20 W Hz" 1 . The radio power of RS2 
is comparable with those of the VLBI sources as- 
sociated with the nucleus A in Arp 299 while the 
power of RSI is nearly a factor ten higher, which 
makes it 1250 times the 8.4 GHz radio power o f 
Cas A of 6 x 10 17 W Hz^ 1 ((Colbert et al.lll994h . 
Both RSI and RS2 display a radio light curve with 
a long-term rise (Figure [5]). The upper limit for 
the source size of RSI at the highest frequency 8.4 
GHz is 1.5 x 0.6 pc (Table [5]) and the source is 
not sufficiently resolved by our EVN synthesized 
beam. Also, different EVN beams resulting from 
different (u,v) coverage at each epoch make it diffi- 
cult to make a reliable comparison of the intrinsic 
source sizes. 

Speculation that RSI could be the true south- 
ern nucleus, instead of Nl, may be ruled out 
because RSI was not detected in the 8.4 GHz 
VLBA observations (1999-2001) and only ap- 
peared in 2003. Also its location is significantly 
offset from the positions of Nl in the 2.4 and 
1.7 GHz VLBA data. Similarly, the possibility 
that RSI represents ejecta from the active nu- 
cleus Nl may be ruled out because its offset (11 
pc) from Nl remained unchanged during 6 years 
and because RSI is far from the known outflow 
structures seen in HI and OH southwest of Nl 
( Baan. Hagiwara fc Hofnerl(2007h . The positional 
difference of 10 mas (re 50 pc) (see Table [2]) be- 
tween RS2 at 1.6 and 5.0 GHz and the source SI at 
1.7 and 2.3 GHz in VLBA observations at epochs 
in 2000 and 2001 suggests that they are also dif- 
ferent sources. 

The radio sources RSI and RS2 have remained 
detectable for more than 6 years since 2003. An 
important question is whether these sources are 
supernovae (SNe, interacting with the circumstel- 
lar matter) or supernova remnants (SNRs, inter- 
acting with the dense ISM). There is no evidence 
for expanding supernovae shells with the highest 
beam size of ~6 mas (~3 pc) both in RSI and RS2 
between the two epochs from our data (see Ta- 
bled]). The range of radio luminosity of if) 20 - 8-21 - 6 
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W Hz 1 of RSI and RS2 is similar to the radio lu- 
minosity obtained for observe d Type Ib/c or Typ e 



2002) 



II radio supernovae (RSNe) (|Weiler et al 
Thus, it is possible that RSI and RS2 are SNe 
or SN-SNR transition objects. Our data is insuffi- 
cient to distinguish between the two cases and fur- 
ther sensitive VLBI monitoring would be required. 
The VLBA data (1999-2001) also displays two 
similar but weaker RSN or SNR candidates RS3 
and RS4, the latter of which has a spectral index 
of a. 1.7-2. 4 ~ 0.35. 

Many nearby galaxies, such as the LIRG M 
82, and ULIRGs Arp 220 and Arp 299, display 
compact radio sources that are evide nce of ongo- 
ing ( ci rcum-)nuclear st a rform a tion ([Smith et al 
1998 : iMcDonald et all [200ll iNeff et all 12004 ; 



at milliarcsecond resolution, that are consistent 
with th e earlier VLBI images obtai ned with the 



Parra et al.l l2007t IConwav et all 120071: lUlvestad 



2009). VLBI observations of the merging galaxy 



Arp 220 have detected a total of 18 compact radio 
sources within the western nucleus of the galaxy, 
over half of which have radio properties that are 
consistent with Type-II supe rnovae interactin g 
with the surrounding medium (jParra et al. 2007 ) . 
Likewise, the nuclei in Arp 299 show 30 compact 
radio sources, 25 of which are associated with the 
northeastern nucleus A (with an AGN candidate) 



and are spread over a region of 30 pc ([Neff et al 
20041 : Illlvestadll2009h . 



A comparison of the SN— SNR sources in 
NGC 6240 with those found in other sources shows 
that the radio powers of RSI and RS2 are equiva- 
lent to tho se of the most powe rful sources found in 
M82 (see iFenech et al.l l2008h . although they are 
less powerful than those in Arp 220. As a result the 
upper limits for the sizes of RSI and RS2 lie in the 
upper range of the general relation between radio 
luminosity and diameter obser ved for Galactic and 
extragalactic SNR sourc es (jBerezhko and Volk 
120041 (Batejat et al.l EhoI ) . This would suggests 
that the environmental conditions in the LIRG 
NGC 6240, and possibly the star formation initial 
mass function, cannot yet be distinguished from 
those of the most luminous FIR galaxies. 

4. Summary 

We have conducted multi-frequency EVN ob- 
servations of the nuclear region of the merging 
galaxy NGC 6240 at two epochs. The new VLBI 
maps reveal the double radio nuclei of NGC 6240 



VLBA (jGallimore fc Beswick 1 120041 ) . The radio 



spectra from both nuclei suggest a spectral turn- 
over between 2 and 3 GHz. The spectrum of the 
southern nucleus Nl may be explained in terms of 
free-free absorption, although this explanation is 
still limited by having source flux densities with in- 
sufficient frequency coverage at each epoch. There 
is no clear interpretation for the spectrum of the 
northern nucleus N2. 

Questions still remain about the true nature of 
the two radio nuclei and whether they both con- 
tain a radio-quiet AGN, a simple starburst, or 
a composite with an AGN and a circum-nuclear 
starburst. We suggest that the southern nucleus 
hosts an AGN and a circum-nuclear starburst, as 
evidenced by the X-ray data and the radio sources. 
The association of the H2O maser with the nuclear 
source is still unknown. However, it is not clear 
that our data, together with earlier VLBA mea- 
surements, confirm the presence of an AGN at the 
northern nucleus. 

The radio components RSI and RS2 could be 
interpreted as long-lived radio supernovae that re- 
sult from ongoing circum-nuclear star formation 
at Nl. The radio spectrum of both sources is rel- 
atively flat and their location remains unchanged 
within error over about 6 years. More radio su- 
pernovae may have been detected in earlier ob- 
servations around the active nucleus of Nl, which 
groups NGC 6240 together with other well-known 
starburst nuclei that display RSNe or SNRs in 
their nuclei such as M 81, M82, Arp 220 and Arp 
299 (|Bartelll2009h . Radio interferometric obser- 



vations are a powerful method to identify ongoing 
nuclear starformation and detect radio supernovae 
and supernova remnants in extragalactic sources. 
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Table 1: EVN Observations of NGC 6240 



Epoch 


Observing 


Frequency 


Telescopes 


Data rate 


Duration 


On-source 




Date 


(GHz) 




(Mbit s" 1 ) 


(hrs) 


(hrs) 


1 


30 Oct 2003 


4.97 


Ef,Hh,Mc,Nt,On,Wb 


256 


7.9 


4.2 


1 


10 Nov 2003 


1.64 


Ef,Hh,Jb,Mc,Nt,On,Wb 


256 


8.2 


3.6 


2 


16 Jun 2009 


4.96 


Ef,Jb,Mc,Nt,On,Sh,Tr,Ur,Wb 


512 


8.9 


5.7 


2 


17 Jun 2009 


2.24 


Ef,Mc,Nt,On,Sh,Ur,Wb 


512 a 


8.7 


5.8 


2 


17 Jun 2009 


8.36 


Ef,Mc,Nt,On,Sh,Ur 


512 a 


8.7 


5.8 



"Total recording rate at both S and X band. 



Table 2: 


VLBI Imagin 


g Parameters 






Epoch 


Frequency 


Rms (Uniform-weight) 


Beam size (FWHM) 


PA 




(GHz) 


(mjy beam -1 ) 


(mas) 


(°) 


1 


1.6 


0.13 


30.8 x 17.8 


50.9 


1 


5.0 


0.083 


10.5 x 4.58 


50.1 


2 


2.2 


0.55 


10.97 x 7.48 


70.1 


2 


5.0 


0.055 


6.11 x 4.93 


72.1 


2 


8.4 


0.11 


3.42 x 2.77 


70.0 



"The rms noise levels near the reference position used in the JIVE VLBI correlation. 



Table 3: Positions of VLBI Components in NGC 6240 



Component 


Coordinate 


Accuracy 




a(J2000) 
16 h 52 m 


<5(J2000) 


Aacos<5 T 




Frequency 


02° 24' 


(mas) 


(mas) 


Northern nucleus (N2) 


S 








1.6 GHz 


58.9242 


04.779 


0.9 


0.7 


2.2 GHz 


58.9240 


04.777 


1.2 


0.9 


5.0 GHz (Oct 2003) 


58.9240 


04.776 


0.6 


0.6 


5.0 GHz (Jun 2009) 


58.9240 


04.776 


0.6 


0.6 


8.4 GHz 


58.9240 


04.776 


0.6 


0.6 


Southern nucleus (Nl) 


1.6 GHz 


58.8903 


03.351 


1.1 


1.1 


5.0 GHz (Oct 2003) 


58.8901 


03.348 


0.7 


0.6 


5.0 GHz (Jun 2009) 


58.8902 


03.350 


0.6 


0.6 


Radio source 1 (RSI) 


5.0 GHz (Oct 2003) 


58.8886 


03.326 


0.7 


0.6 


5.0 GHz (Jun 2009) 


58.8885 


03.326 


0.6 


0.6 


8.4 GHz 


58.8885 


03.327 


0.6 


0.6 


Radio source 2 (RS2) 


1.6 GHz 


58.8991 


03.607 


2.0 


1.3 


5.0 GHz (Oct 2003) 


58.8993 


03.602 


0.8 


1.3 


5.0 GHz (Jun 2009) 


58.8994 


03.602 


0.6 


0.6 



tThe position error in right ascension (Aacos5) and declination (A5). 
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Table 4 



Component Flux Densities of NGC 6240 





Epoch 1 (EVN;Oct- 


-Nov 2003) 






Epoch 2 (EVN;June 2009) 




MERLIN 


Name 


Sl.6 


S5.0 


«i:8 




S2.2 


S5.0 


Ss.4 




S5.0 


(i) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


Nl 


2.09 


2.63 


0.20 




< 2.75 


2.85 


<0.55 


<-3.2 


36.8 


N2 


3.95 


5.99 


0.37 




4.59 


5.79 


1.78 


-2.3 


16.1 


RSI 




2.34 






<2.75 


3.00 


3.57 


0.34 




RS2 


1.07 


0.38 


-0.91 




<2.75 


0.59 


<0.55 


<-0.14 




VLBA 




Epoch 1 (Aug 1999) 




Epoch 2 (Apr 


-May 2000) 




Epoch 3 (Dec 2000-Jan 2001) 






Si. 7 


S2.4 




Si. 7 


S2.4 




Si. 7 


S2.4 




(10) 


(11) 


(12) 




(13) 


(14) 




(15) 


(16) 




SI 




0.63±0.13 




0.62±0.06 


0.59±0.10 




0.40±0.08 


0.47±0.08 




RS3 
















0.46 




RS4 








0.44 


0.29 











Note. — Flux densities are in mjy. Col. (1): Source names from Baan, Hagiwara, Sz Hofncr (2007) for Nl (southern nucleus) 
and N2 (northern nucleus). RSI and RS2 arc newly detected sources in these EVN observations. Flux density errors measured 
by the EVN arc estimated at ~ 10%. Col. (2), (3): The 1.6 GHz and 5.0 GHz flux densities at Epoch 1 as measured by the 
EVN. Col. (5)-(7): The 2.2 GHz, 5.0 GHz, and 8.4 GHz flux densities at Epoch 2 as measured by the EVN. Col. (4), (8): The 
spectral indices between 1.6 GHz an d 5.0 GHz (Ep o ch 1) and between 5.0 GHz and 8.4 GHz (Epoch 2). Col. (9): Data taken from 
MERLIN obser vations at 5.0 GHz in l Beswick ct al. (2001). All flux density upper limits arc 5 a. Col. (10): Data taken from VLBA 
observations in Gallimorc & Bcswick (2004). RS3 (~0 // .02 northwest of the southern nucleus Nl) and RS4 (^0 // .03 northwest of 
the norther n nucleus N2) are newly defi ned in this article. Col. (11)-(16): The 1.7 GHz and 2.4 GHz flux densities as measured by 
the VLBA dGallimore fc Beswick~ll2004|) . 
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Table 5 



Source sizes of Nl, N2, RSI, and RS2 







Nl 






N2 






RSI 






RS2 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Frequency (Epoch) 


(mas 2 ) 


(pc 2 ) 


(°) 


(mas 2 ) 


(pc 2 ) 


(°) 


(mas 2 ) 


(pc 2 ) 


O 


(mas 2 ) 


(pc 2 ) 


(°) 


1.6 GHz (Oct 2003) 


17.6x12.6 


8.4x6.0 


43 


39.3x13.5 


18.7x6.4 


5 








15.8x7.0 


7.5x3.3 


169 


5.0 GHz (Oct 2003) 


9.2x5.9 


4.4x2.8 


137 


19.5x6.6 


9.3x3.2 


98 


9.2x3.9 


4.4x1.9 


80 


9.2x6.1 


4.4x2.9 


140 


5.0 GHz (Jun 2009) 


7.6x5.2 


3.6x2.5 


161 


12.4x3.4 


5.9x1.6 


101 


3.6x1.9 


1.7x0.9 


70 


2.5x2.9 


1.2x1.4 


128 


8.4 GHz (Jun 2009) 








6.6x1.4 


3.1x0.7 


85 


3.2x1.3 


1.5x0.6 


58 









Note. — Source sizes are upper limit values and they are deconvolved using the VLBI synthesized CLEAN beam. Col. (1): The observing frequency and 
Epoch. Col. (2), (5), (8), (11): Major axis X minor axis values, in milli-arcsecond (mas). Col. (3), (6), (9), (12): Major axis X minor axis values, in parsccs (pc). 
Col. (4), (7), (10), (13): The position angles of the sources. 



Table 6: Properties of the VLBI sources Nl, N2, RSI, and RS2 


Source 


Tb,5GHz 


P (5.0 GHz) 


Tb, 8.4GHz 


P (8.4 GHz) 




(x 10 6 K) 


(x 10 21 W Hz" 1 ) 


(x 10 6 K) 


(x 10 21 W Hz" 1 ) 


Northern nucleus (N2) 


10.0 


2.9 


4.9 


1.7 


Southern nucleus (Nl) 


5.2 


1.7 


<1.5 


<0.70+ 


Radio source 1 (RSI) 


31.6 


3.9 


21.8 


4.1 


Radio source 2 (RS2) 


5.9 


0.75 


<1.5 


<0.70t 



Note. — These values are calculated from data obtained in June 2009. f5 a upper limit 



Table 7 



Modeling fit parameters for the nuclei Nl and N2 



Parameter 


Value 


Free-free absorption model 


Nl 


N2 


So 


25.4 ± 35.7 mjy 


30.4 ± 60.8 mjy 


a 


-1.4 ± 0.8 


-1.1 ± -1.1 


Tff 


5.3 ± 3.5 


4.5 ± 4.9 


Reduced \ 2 


0.59 


3.8 


Synchrotron self-absorption model 


Nl 


N2 


So 


0.39± 0.07 mjy 


0.74± 0.22 mjy 


a 


-3.2 ± 2.3 


-2.40 ± 1.97 


T B 


1349 ± 5110 


400 ± 1326 


Reduced x 2 


0.43 


2.9 
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02 24 05.0 — 



O 

5 



o 

LU 
Q 



04.4 — 



04.0 



03.8 



03.4 



03.2 — 



03.0 




16 52 58.95 58.94 58.93 58.92 58.91 58.90 58.89 58.88 58.87 58.86 
RIGHT ASCENSION (J2000) 



ASCENSION (J2000) 



Fig. 1. — The radio struct ure of the nuclear region of NGC 6240 at 5 GHz as obtained by MERLIN 
( Gallimore fc Beswick 20041 ) with the location of compact components detected by the new 5 GHz EVN 
observations in October 2003 (Epoch 1) as presented in the right-hand-side frames. The synthesized beam 
(FWHM) of the observations is plotted in the left bottom corner of all maps. The two radio nuclei Nl 
(southwest) and N2 (northeast) have been labelled fol l owing the designation of Colbert et al. (1994) and 
Baan et al. (2007). Note that in lGallimore fc Beswick I (|2004l ) the northern nucleus is labeled as Nl and the 
southern nucleus labeled as S. The cross indicates the positions of radio source RS2 that appeared in our 



observation. The contours in the MERLIN image are at 0.5 mJy beam 
of the three VLBI components are given with Figures 2, 3, and 4. 



-1, 1, 2, 4,..., 32. The contours 
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RIGHT ASCENSION (J2000) 
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Levs = 1.10DE-04 " (-4, 4, 10, 16, 22) 



Fig. 2. — Radio continuum images of the southwestern nucleus Nl of NGC6240 and the compact radio 
source RSI, which is located southwest of Nl. The synthesized beams are plotted at the bottom-left corner 
on each panel, (a) The 1.6 GHz image of Nl observed in Epoch 1 in 2003. (b) The 5.0 GHz images of Nl 
and RSI observed in Epoch 1 in 2003. (c) The 5.0 GHz images of Nl and RSI observed in Epoch 2 in 2009. 
(d) The 8.4 GHz image of RSI obtained in Epoch 2. 
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Fig. 3.— Radio continuum images of the northeastern nucleus N2 of NGC 6240. The synthesized beams are 
plotted at the bottom-left corner on each panel, (a) The 1.6 GHz image observed in Epoch 1 in 2003. (b) 
The 2.2 GHz image observed in Epoch 2 in 2009. (c) The 5.0 GHz image observed in Epoch 1 in 2003. (d) 
The 5.0 GHz image observed in Epoch 2 in 2009. (e) The 8.4 GHz image obtained in Epoch 2 in 2009. 
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Fig. 4. — Radio continuum images of the radio source RS2 that is located northeast of the southern nucleus 
Nl as indicated in the lower angular resolution map obtained by the 5 GHz MERLIN (Figure [T]). The 
synthesized beams are plotted at the bottom- left corner on each panel, (a) The 1.6 GHz image obtained in 
Epoch 1 in 2003. (b) The 5.0 GHz image of Epoch 1 in 2003. (c) The 5.0 GHz image of Epoch 2 in 2009. 
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Fig. 5. — Radio spectra of the two nuclei of NGC 6240. 

Left frame: The spectrum of the northern nucleus N2. The solid lines represent the free- free absorption 
(FFA) model as fitted to the Epoch 2 data and 1.6 GHz data in Epoch 1, assuming no significant flux 
density variation of N2 between the two observing epochs. The dotted line represents the synchrotron self- 
absorption (SSA) model as fitted to the same data as the case of FFA. 

Right frame: The spectrum of the southern nucleus Nl. The solid and dotted lines represent the FFA model 
and SSA model as fitted to data of the two EVN epochs including upper limit values, assuming no significant 
flux density variation of Nl during these two epochs. 

The data points: The filled circles at 2.2, 5.0, and 8.4 GHz are data points from the EVN during Epoch 2 
(2009). The open circles at 1.6 and 5.0 GHz are data points from the EVN during Epoch 1 (2003). The 
open squares at 1.7 and 2.4 GHz are the data points from the VLBA at Epoch 2 (April to May 2000), and 
filled squares at 1.6 and 2.4 GHz are from the VLBA data of Epoch 3 (December 2000 to January 2001). 
The triangles denote the 5 a uppe r limit values obtained at 2.2 and 8.4 GHz with the VLBA in August 1999 
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Fig. 6. — Radio light curves of RSI, RS2, and SI. Filled circles denote the 5.0 GHz flux density of RSI, and 
an open circle denotes the 8.4 GHz flux density of RSI. A filled square shows the 1.6 GHz flux of RS2, and 
the open squares denote the 5.0 GHz flux of RS2. Filled and open diamonds denote the 1.7 GHz and 2.4 
GHz flux of SI from the earlier VLBA data. A open triangle indicates the upper limit (5 a) of the 8.4 GHz 
earlier VLBA data. 
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